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We model electric field noise from fluctuating patch potentials on conducting surfaces by taking 
into account the finite geometry of the ion trap electrodes to gain insight into the origin of anomalous 
heating in ion traps. The scaling of anomalous heating rates with surface distance, d, is obtained for 
several generic geometries of relevance to current ion trap designs, ranging from planar to spheroidal 
electrodes. The influence of patch size is studied both by solving Laplace's equation in terms of the 
appropriate Green's function as well as through an eigenfunction expansion. Scaling with surface 
distance is found to be highly dependent on the choice of geometry and the relative scale between the 
spatial extent of the electrode, the ion-electrode distance, and the patch size. Our model generally 
supports the d"** dependence currently found by most experiments and models, but also predicts 
geometry-driven deviations from this trend. 
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I. INTRODUCTION 

Metallic conductors are widely considered to be ideal 
electric equipotentials. However, from the perspective 
of precision measurements and quantum information sci- 
ence, real metals are surprisingly poor equipotentials. 
Not only do metallic surfaces exhibit static, inhomoge- 
neous, microscopic potential patches [IH3]j but even 
more interestingly, the electric fields produced by real 
conductors exhibit time- varying fluctuations. These ef- 
fects have profound implications for many modern exper- 
iments: the static patches produce static electric fields, 
which influence precision measurements ^E], while the 
time-varying fluctuations, or electric field noise, limit 
progress on single spin detection [7], nanomechanics [S]- 
[TO] , detection of weak forces [SJ [TTl [T^] , and quantum in- 
formation processing with trapped ions |13H16| . As such, 
many fields of inquiry would be advanced if the physical 
origin of these phenomena were to be understood. 

A few physical models for patch potentials and elec- 
tric field noise exist. The static component of patch po- 
tentials has been studied in great detail and is believed 
to be related to work function differences between crys- 
tal facets and to surface adsorbates [TJ [3] . Electric field 
noise, however, is less understood, with two main candi- 
date models. One model explains that this noise is a con- 
sequence of the finite electrical resistance of the metallic 
electrodes and associated external circuitry used in ex- 
periments, otherwise known as Johnson noise jl3lll7fE0] . 
A second model attributes the noise to patch potentials 
with a fluctuating component [131 121j , and is motivated 
by experimental evidence which substantiates a 1/w fre- 
quency scaling [T31 [TH [TB] , with an origin rooted in ther- 
mally activated processes [3 [TH (THl [22] . The fluctuating 
patch potential noise model is the focus of this paper. 

Differentiating between noise models requires accurate 
and precise measurements of noise. Depending on the 
probing frequency and the distance from the surface at 



which the electric field noise is measured, very different 
methods may be used, such as cantilevers [BHTO]. Kelvin 
probes [I1II2], or trapped ions [I3HI6] . In particular, 
atomic ions, cooled to the motional ground state of their 
trapping potential, exhibit remarkable sensitivity to elec- 
tric field noise. Experiments can discern changes in the 
motion of the ions to within a fraction of single quanta 
due to extremely narrow electronic transitions in ions, 
which can feature fractional linewidths of about one part 
in 10^^. This sensitivity has revealed an anomaly in the 
electric field produced by conductors: it is significantly 
noisier than predicted by Johnson noise, leading to a phe- 
nomenon known as anomalous ion heating |13H16l I22j . 

Anomalous heating in ion traps has been studied, e.g., 
by Turchette et al. 13] by modeling the ion as a single 
charge at the center of a conducting sphere [Fig. [l] (a)], 
which is a reasonable approximation to, e.g., the hy- 
perbolic radio frequency (RF) Paul trap [Fig. [I] (b)]. 
Their work showed that the power of electric field noise 
should scale with ion-electrode (conductor) distance d as 
dr" with a scaling exponent a = 2 for Johnson noise, 
and a = 4 for patch potential noise. From collective 
data gathered by several independent ion trap experi- 
ments [15], a pattern in support of the patch potential 
model has emerged; however, comparison between exper- 
iments is fraught with uncertainty as anomalous heating 
is known to be strongly influenced by trap preparation 
[T3l [22] . A systematic study of the distance scaling in a 
single ion trap was performed by Deslauriers et al. |14j . 
who found a = 3.5 ± 0.1. While in reasonable agree- 
ment with the model of fluctuating potential patches, 
the needle-shaped electrode geometry of their ion trap 
[Fig. [T] (c)] bore little resemblance with that of the ion- 
in-a-sphere model of Turchette et al. 

It is thus of interest to model the scaling of anoma- 
lous heating in the patch potential model for different 
geometries. In this respect, the Turchette et al. model is 
complimented by the model of Dubessy et al. [21], who 



2 




FIG. 1. Examples of electrode geometries relevant to ion trap- 
ping, d is the nearest ion-electrode surface distance. Elec- 
trodes are at either RF potential (shaded) or DC potential 
(unshaded), a) lon-in-a-sphere model used in Ref. [2]. b) 
Cross sectional view of the hyperbolic Paul trap, c) Needle 
electrode Paul trap, d) Surface electrode Paul trap. 



considered an infinite planar surface. The infinite planar 
geometry is representative of the technologically impor- 
tant surface electrode ion trap [53H2S], where the ion is 
trapped above the surface of two-dimensional electrodes 
[Fig. [I] (d)], as well as of the situation encountered in 
non-contact friction measurements with cantilevers [S]. 
Dubessy et al. |21j recognized the importance of mod- 
eling finite correlations between potential patches and 
introduced a geometric factor, for the patch correla- 
tion length, which led to a prediction of a strong de- 
pendence of a on this length scale. Specifically, they 
obtained a = 4 for c? ^ ^ and a = 1 for o? <C Ci which 
allowed for an impressive connection between data from 
ion trap to cantilever experiments across more than three 
orders of magnitude in scale. 

Here, we study the influence of electrode geometry on 
patch potential driven electric field noise in ion traps by 
extending previous theoretical work [T31 HH US] , to allow 
for a consideration of an arbitrary geometry and patch 
correlation function. By applying the method of eigen- 
function expansions, we are led to an intuitive under- 
standing of the effects of arbitrary patch distributions 
and we apply our model to explicit examples of simple 
finite geometries where the scaling exponent a is eval- 
uated. In these simple finite geometries, we emphasize 
the limiting cases of infinite and infinitesimal patch sizes. 
This enables us to establish bounds on a across a range 
of relative scales of patch size, ion-electrode distance, and 
electrode dimensions. 

The specific electrode geometries we consider include 
the infinite plane as well as finite spheroidal electrodes. 
The planar geometry is representative of typical surface 
electrode ion traps and is presented, in part, 

to connect with previous works [HIISS]. Spheroids, on 
the other hand, are examples of finite geometries, where 
limiting cases can be chosen to approximate both finite 
planar electrodes and the needle trap of Deslauriers et 
al. [H]. 

We find that a assumes a wide range of values de- 



pending on the relative scales of patch size, ion-electrode 
distance, and spatial extent of the electrode. Near elec- 
trode surfaces, the scaling exponent is bounded within 
the range < a < 4, depending on patch size. In typical 
ion trap configurations with a large spatial extent of elec- 
trodes relative to an ion height of d ^ 100 /im and with 
small patch sizes of dimension ^ 1 /im, a approaches 4, 
in agreement with the models of Turchette et al. and 
Dubessy et al. As the distance, d, is decreased we find 
that a decreases, as was also determined by Dubessy et 
al. For the special case of the planar hole trap |171 [55], 
with the ion suspended at the center of a hole in the 
RF electrode, we find that the bounds are narrowed to 
2 < a < 4 but the a = 4 scaling is retained in the limit of 
typical trap dimensions. In the case of spheroidal shapes, 
a converges to either 4 or 6, when the ion is far from 
the spheroid, depending on whether the mode of motion 
considered is normal or transverse to its surface. For the 
limiting case of a needle in the intermediate regime of 
d ~ a, geometry and patch size dependence may allow 
for scaling consistent with the a = 3.5 value for the nor- 
mal mode of motion, as found by the Deslauriers et al. 
experiment [14 . 

Our formalism and its results are presented in the fol- 
lowing. We begin, in section [ll] by briefly reviewing the 
general theory for heating rates in ion traps due to fluc- 
tuating patch potentials in the framework of Laplace's 
equation and its Green's function solution. In section [ITT] 
we study the influence of electrode geometry on the scal- 
ing of heating rates for various generic geometries of rel- 
evance to current ion trap designs. Finally, in section [IV| 
we summarize key results of our work and propose direc- 
tions for further studies, where the analytical methods 
developed here may aid in numerical simulations of non- 
generic geometries to model more realistic experimental 
scenarios. 



II. MODEL 



We analyze the scaling laws behind heating rates as 



follows: In section II A the dependence of single ion heat- 
ing rate on fluctuating electrical flelds is determined. In 
section |II B| it is assumed that these fluctuating electric 
flelds originate from potentials on some conducting sur- 
face with a certain geometry. In section II C the patch 



potential approximation is implemented, which allows all 
information about electrode and patch geometry to be 
described by a single geometric factor, which we shall 
denote A. In section II D assumptions about patch sizes 



are introduced to evaluate upper and lower bounds for 
the scaling laws. Finally, in section |II E| the scaling ex- 
ponent a is derived from the geometric factor. 
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A. Ion heating rate 



Our treatment of electric field noise in ion traps as- 
sumes that the ion is confined in a harmonic potential in 
three dimensions, where the strength of the potential is 
quantified by frequencies Wfc, with k denoting the prin- 
ciple axes of the potential. Such confinement may by 
achieved either by a combination of static electric and 
magnetic fields, as in the case of the Penning trap, or by 
a combination of static and time-varying electric fields, 
as in the case of the Paul trap ^29| . Our focus is on the 
Paul trap, of which a few examples are shown in Fig. [T] 
b)-d), but our formalism is quite general, and applies to 
any scenario where a particle of mass m and charge q is 
held in a harmonic potential near a conducting surface. 
For this reason, we shall not review the subject of ion 
confinement here; we refer the reader to Refs. [211 [SSI US] 
for detailed treatment of that subject. 

We assume the ion is initially is the quantum mechan- 
ical ground state of the harmonic oscillator potential and 
consider the effect of a fluctuating electric field with com- 
ponents Ek{t) at the ion location. The heating rate is 
now defined as the rate at which this field induces tran- 
sitions from the ground state to the first excited state 
and can be evaluated via first order perturbation theory 
to [131 



Amhiij] 



(1) 



Here 
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Ssd^k) = 2 / dTe''^''^Ek{t)Ek{t + r)) (2) 



is the power spectrum of the electric field noise and (• • • ) 
denotes time averaging. Cross coupling between the noise 
and the RF drive field at frequency 17 occurs in Paul 
traps, but has negligible effect as long as (wfc/fi)^ <C 
1 1131 and is hence omitted from our treatment. 



B. Electric field noise spectral density 



Literature on patch potentials [TH2| suggests various 
origins such as crystal grain boundaries, which give rise to 
variations in local work functions on the conductor, and 
adsorbed elements on the surface, which may alter the 
surface potential locally as well. Sub-surface defects have 
also been considered; however, measurements with ions 
trapped above superconducting surfaces that can provide 
shielding against fields from such defects, suggests that 
their influence is weak |30j . In the following, we thus 
focus on a description based on a surface effect. 

The fluctuating electric field Ek{t) (Eq. [2]) can be 
obtained by solving Laplace's equation with Dirichlet 
boundary conditions specified by the potential </>, which 
is determined by fluctuating patch potentials on a sur- 
face (7 of interest. In the Green's function formalism for 
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FIG. 2. Graphic illustration of Eq. |3] for a given, arbi- 
trary conducting surface cr, here represented by the shaded 
cylinder. The electrostatic potential $(r) at point r is de- 
termined by adding the convolution of some charge density 
p(r') with the Green's function G(r,r') and the convolution 
of the surface potential ^(r") with the surface Green's func- 
tion -a„»G(r,r"). 



a time-dependent electrostatic charge density, p(r, t), the 
potential is given by 

$(r, t) = G(r, r')p(r', t)dv' - ^ 0(r", ^^^rfr", 

(3) 

where G(r,r') is the Green's function of Laplace's equa- 
tion for some arbitrary conducting surface. V is the vol- 
ume of integration, and n" is normal to the surface of 
integration a at some point r" G a (Fig. [2]) . We have ap- 
plied here the quasi-static approximation as typical ion- 
electrode distances of d ~ 100 /im are much less than the 
relevant wavelength of electric fields at the secular mo- 
tion frequencies, which are typically of order ~ 10 MHz. 

As we are most interested in surface effects, p repre- 
sents charges adsorbed onto the surface, whereas (/) rep- 
resents patch potentials. However, the electric potential 
of a charge very close to a conducting surface is indistin- 
guishable from an appropriately chosen patch potential 
as long as the ion-surface distance is much larger than 
the adsorbed charge-surface distance. Hence, for sim- 
plicity, we set p = and treat all sources in terms of 
effective patch potentials. The equivalence between sur- 
face and free-space sources implies that, for sources </) on 
the surface, — 9„'G(r, r')|r'Go- is analogous to G(r,r') for 
sources p in free space. We therefore define the surface 
Green's function Gc{v,y') = — 9„'G(r, r')|r'e<j. In this 
notation, the temporal correlation function of the elec- 
tric field along the /c*'^ mode is 

{Ek{t)Ek{t + t)) - (Vfc$(r,i)Vfc$(r,t + r)) 

(</.(r',t)0(r",i + r)) 



X (VfcG<,(r,r')-VfcG,(r,r") dr'dr" 



(4) 



where Vfc is the k component of the gradient operator. 
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Inserting into Eq. [2] we find 

SEju;k)^2 I f F[(0(r',t)0(r",t + T))] 

cr" ^ a' 

X fVfeG,(r,r')VfcG,(r,r")Vr'dr", (5) 



where F [• • • ] = [• • • ] e^^'^dr denotes the temporal 
Fourier transform. 



C. Patch potential approximations 

In order to evaluate Eq. [5j some assumptions need to 
be made about the nature of the sources. One common 
assumption, such as in Ref . |21) , is that the temporal and 
spatial variation of the sources decouple. In this approx- 
imation, the sources are described by a superposition of 
separate patches, where the i^^ patch is described by 
a time dependent function Vi{t) and an effective spatial 
extent Xi{^) ~ hence the term patch potential. Thus, 



N 



(r,t)=^F,(i)x.(r), 



(6) 



1=1 



where the normalization of Xr{'^) is with respect to the 
area A of the surface a: 



N 



(7) 



A measure of the average patch size may be obtained 
through the spatial correlation function 



N 

C(r',r")=^X.(r')-X.(r"), 

4=1 



with normalization given by 

/ C(r,r')drdr' 

J a 



N ' 



(8) 



(9) 



This equality is satisfied for both patches that overlap 
maximally and disjoint patches with zero overlap. The 
former case is unphysical for naturally occurring patches, 
but could be reached experimentally by applying strong 
electric noise from a single source across all surfaces. 
Lastly, the sources are assumed to share the same spec- 
tral distribution, R{lo), but are also uncorrelated with 
each other, such that 



F[{V,{t)-Vjit + T))] = RiLu)S,. 



(10) 



The noise spectrum has been observed by several 

experiments to resemble a l/w scaling [I3UI6] . although 
theoretical models exists which supports a frequency de- 
pendence extending to both l/co^-^ [3T] and l/w^ [25] . 
depending on the assumptions used. 



Given our assumptions above, the frequency depen- 
dence in the expression for the spectral density of Eq. [5] 
is conveniently separated from the spatial dependance 
and becomes 



N 



.i=l 

X (VfcG,(r, r') VfeG, (r, r")) dr'dr" 
= 2i?K)Afc(r), 
where we have defined the geometric factor 

A,(r)= / / G(r',r") 

X (VfcG,(r, r')VfeG,(r, r")) dr'dr". 



(11) 



(12) 



which encompasses information about patch sizes 
through G(r', r") and the surface geometry through the 
surface Green's function Go.(r,r'). 



D. Patch size 



It is evident from Eq. 12 that Ga(r, r') and G(r',r") 



are generally not separable in the expression for the 
power spectrum of the electric field noise. Furthermore, 
analytic expressions for the geometric factor do not exist 
for most choices of geometries or correlation functions. 
We deal with this issue in two ways: One, by consider- 
ing limiting cases of the correlation function, and two, by 
considering an eigenfunction expansion of the correlation 
function. 

The first approach takes two limiting cases of patch 
size relative to the spatial extent of the surface where 
the problem simplifies, which provides some useful in- 
tuition about the behavior encountered in intermediate 
regimes. The limiting cases are the Infinite Patch (IP), 
which corresponds to the trivial case of a conductor held 
at uniform potential, and the Point Patch (PP), repre- 
sented by an infinitesimally small patch disjoint from all 
other patches. Formally, 



G(r',r") 



1 

7V" ^ 



IP, 
PP. 



(13) 



The choice of AjN as the coefficient for the PP limit 
satisfies the normalization in Eq. |9j and N is formally 
given by A ■ limr-j.o (5^(r). In these limits, the geometric 
factor takes the following forms: 

A,(r)-/I^'^'^^^^(^'^')^^'l' (14) 

The physical intuition behind Eq.[l4]is that within a large 
(infinite) patch, noise from all sub-elements of the patch 
are correlated and the fields add coherently, while for 
small (point) patches, the contributions from individual 
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C(r,r') Depiction of patches 



coefficients Ci, in the IP and PP Hmits are 



5{r - r') 




IP 



PP 



FIG. 3. Patch sizes and associated correlation functions, rang- 
ing from the Infinite Patch (IP) to the Point Patch (PP) from 
top to bottom. 



points are uncorrelated and the fields are added incoher- 
ently by summing the intensities. As long as patch sizes 
are considerably smaller than the variation of the surface 
geometry, this leads us to expect a smaller magnitude 
of noise for smaller patches than for very large patches, 
in agreement with the conclusions of previous models of 
noise originating from patch potentials [13l [21] . 

Our second approach is motivated by the fact that, 
even in the IP and PP limits, it is often not the case 
that the geometric factor can be evaluated analytically. 
Furthermore, it is useful to evaluate the geometric fac- 
tor for some specific correlation function that parame- 
terize patch sizes. To accommodate with these issues, 
we present a formalism for the geometric factor that is 
tractable by numerical methods. We start with the eigen- 
function expansion of the Green's and surface Green's 
functions J32, : 



X, dnfi {r)dr dnf* (r)dr 
^X.a„/.(r)5„/;(r)dr 



IP, 
PP. 



(18) 



The patch correlation functions in Eq. [13] represent two 
opposite extremes; however, Eq. [17] facilitates the study 
of the intermediate regime. Moreover, the latter raises 
the possibility of constructing more complex correlation 
functions out of a relatively small set of basis correlation 
functions that are simultaneously diagonalizable in the 
eigenfunction basis. This approach lends itself well to nu- 
merical modeling and is generally applicable to arbitrary 
geometries and patch correlation functions. By contrast, 
e.g., the choice of C(r,r') cx e"!'-'"'!/'; in Ref. [21] facili- 
tates an analytical solution for an infinite planar surface, 
but not for other geometries. 

Patch sizes can be very dependent on material and sur- 
face preparation; however, values in the range ~ 10 nm to 
^ 1 /im are typically reported by experiments [H El IS] . It 
follows that, at ion-surface distances of c? ~ 100 /im that 
are typical to many ion trap experiments, the exact form 
of the correlation function will not influence the evalua- 
tion of electric field noise. We shall thus not concern our- 
selves with the form of the correlation function, nor with 
the physics governing it, and only use a simple parame- 
terization of intermediate patch sizes by truncating the 
eigenfunction expansion beyond some higher-order term 
in the sum of Eq. |17| Ex plicit examples of this method 
are given in section jlll C and III D 



G(r,r')=^ 



G.(r,r') = -5] 



/.(r)/.r(rO 
A,; ' 



(15) 
(16) 



where /i(r) are eigenfunctions of Laplace's equation with 
eigenvalues Xi and satisfy homogenous boundary condi- 
tions on the surface a. Substituting the result into the 
expression for the geometric factor and using G(r, r') = 
G*(r,r') gives 



Afc(r) = ^c,,Vfc/,(r)V,/;(r), 



/^,/^G(r,r')5„/,(r)a„,/;(r')drdr' 



a,a; 



(17) 



The geometric factor is now expressed as a double sum 
over eigenfunctions with expansion coefficients Cij , which 
represent projections of the correlation function onto 
eigenfunctions of the surface geometry. The average 
patch size can then be modeled by specifying a correla- 
tion function that parameterizes the patch size Q (Fig.[3]), 
and projecting onto these eigenfunctions. The expansion 



E. Scaling with surface distance 

The focus of this work is on the scaling of the power 
spectrum of the electric field noise with surface distance, 
described through the geometric factor Afc(r). Prior work 
on this subject has typically found this to be described by 
some power law with respect to surface distance d, i.e., 
Afe(r) cx . We shall assume such a relationship to 
hold in general, and our parameter of interest, the scaling 
exponent a, is evaluated as 



(19) 



ak{d) = -9in[d] (In Afe(ro + dd) 



where Tq is a reference point on the surface against which 
scaling in some direction d is evaluated. In the following 
we consider various surface geometries and we quantify 
their infiuence on the geometric factor through a in this 
way. Note that if Af^ (r) scales differently for the different 
k modes, it is possible to obtain different values for a 
if, for example, motion normal instead of transverse to 
the electrode surface is considered. While this has been 
considered in past treatments of Johnson noise [20], it 
has largely been overlooked by patch potential models. 



FIG. 4. Coordinate system used in the Green's function for 
the infinite plane in Eg. |20[ The black dot represents the ion. 



FIG. 5. Ion (black dot) in a planar hole trap with radius d. 
RF is applied to the entire plane. 



III. ELECTRODE GEOMETRY DEPENDENCE 

We now proceed to consider the electrode geometry of 
the ion trap. The focus is on generic structures, including 
finite geometries, often found in experiments. We will 
consider planar electrodes and special cases thereof with 
hole-in-conductor, as well as spheroidal shapes such as 
spheres and very prolate spheroids as models for needles. 



A. Infinite planar electrode 

We first establish a connection with previous work in 
the field [21] by considering the infinite planar electrode. 
While inherently unphysical, this represents a good ap- 
proximation to many ion trap geometries where the ion 
is confined only ^ 100 /im above a surface extending sev- 
eral mm in both directions of the plane |23j. It is also 
an excellent model for the geometry encountered in non- 
contact friction measurements with cantilevers [B] - a fact 
exploited by Dubessy et al. PT to connect measurements 
of noise in these seemingly disparate systems through a 
single physical model. 

The Green's function for an infinite conductor is ob- 
tained from a straightforward consideration of image 
charges. In cartesian coordinates {x,y,z) where the z 
axis is normal to the surface (Fig. [4|: 



G(r,r') 



1 

47r 



r — r' 



r — r' 



2z'z 



(20) 



When an adsorbate with charge q at r' is brought 
much closer to the surface than the ion-surface dis- 
tance, the corresponding potential seen by the ion is 
indistinguishable from that of a point patch potential 
(j) = qz'S'^{{x' — x)x + {y' — y)y), which is consistent with 
earlier statement that the electric potential of a charge 
very close to a conducting surface is indistinguishable 
from an effective patch potential. The physical origin 
of Vi{t) in this case would then be fluctuations in the 
charge-surface distance z' . 

The infinite (point) patch limit is reached when the 
ion-surface distance is much smaller (larger) than the av- 



erage patch size. By substituting the Green's function of 
Eq. 20 into the expression for Afc(r) in Eq. 12 we obtain 



A^{dz) 



K,y{dz) 







N 327r" 



IP, 

PP, 

IP, 
PP. 



(21) 



Ion trap experiments studying this geometry, such as 
those based on the surface electrode ion trap [Fig. [T] (d)], 
have focused on the point patch limit as ion heights of 
d ^ 100 /im and a trap size of ~ 10 mm, respectively, 
are typical. At these dimensions, patch sizes of ~ 1 /zm 
appear as point patches and the model thus predicts an 
a = 4 scaling. As the ion-surface distance is lowered, one 
expects a to decrease, as was also concluded by Dubessy 
et al. [21 . In principle, a continuous transition from 4 
up to is predicted by our model; however, we note that 
the IP limit is trivial, as an infinite plane at uniform 
potential has no electric field. Furthermore, it is clear 
that, as systems are scaled to small surface distances be- 
low the dimensions of the patches, the exact form of the 
patch correlation function will have a strong impact on 
the geometric factor 



B. Infinite planar electrode with a hole 

This section considers the geometry of the hole trap 
studied previously in e.g. Refs. (171 [55] and depicted in 
Fig. [5] This geometry is also of relevance to large-scale 
quantum information processing where an architecture 
based on an array of such traps has been proposed [33j . 
We model this hole trap as a thin, infinite conducting 
sheet with a circular hole. As the ion is always trapped 
at the center of the hole, the relevant scaling parameter 
becomes the hole radius, which we label here as d. In 
cylindrical coordinates (Fig. [6]) , where 



x — s cos ( 



y ~ s sm 0, z 



(22) 



the Green's function for a thin conductor with a hole 

is m 
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FIG. 6. Coordinates system used for the hole trap. The ion 
(black dot) is suspended at point r. The source is located at 
r'. 



1 ( 1 
1 



, 2 

1 H — arctan — 
TT V r_ 



2 

1 + e— arctan 

TT 



+ 



(23) 



where r, r' lie on the same side of the plane, and where 



2ss' cos ( 



b') + {zTz'y, (24) 



S^ + Z^- d^) + z'2 - d^) 



±Ad^zz' ■ 



Z2 + (s - df Z2 + + rf)2 
^ 1/2 



Z'2 + (s' - d)' Z'2 + (s' + d) 



(25) 



and 



e = sgn [z + z'2 - ^2^ ^ ^, ^_^2 ^ ^2 _ ^2^] ^ (26) 

We place the ion at the center of the hole at r = (0, 0, 0) 
and substitute the Green's function for a conductor with 



a hole (Eg. [23|) into the geometric factor (Eq. 12). Care 
must be taken when evaluating the integral in the geo- 
metric factor as the sharp edge of the hole at s = d results 
in a divergence. To circumvent this issue, one may in- 
troduce a small parameter 5 ^ d and perform the radial 
integration over (d -I- 5, oo). Doing so we obtain 



A,(0) 



A,(0) 




A 1 



;21n^-3))d- 



IP, 
PP. 



(27) 



Note that the logarithmic divergence in d/5 for radial 
scaling does not detract from the main result of a = 4. 

Unlike the infinite plane, the IP limit of does 
not evaluate to zero as only one side of the sheet is held 
at a non-zero potential. It is furthermore interesting to 
note that similar results for the scaling of the geometric 
factor are obtained as for the infinite plane despite the 
significant difference between the two geometries. 



FIG. 7. Approximations to the Deslauriers et al. needle 
trap [14]. Black dot: ion. Shaded electrode: RF. a) Nee- 
dle trap configuration, b) Approximation to a) assuming a 
single needle suffices to describe the scaling of the geometric 
factor with respect to d. c) approximation to b) representing 
the needle by an effective sphere of radius a. 
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FIG. 8. Coordinates system used for the sphere. The ion is 
suspended at point r. The source is located at r' 



C. Spherical electrodes 

We now consider a spherical electrode, which is an ex- 
ample of a finite geometry where placing an ion outside 
the sphere and far from its surface is well defined as well 
as an example of a geometry radically different from the 
two-dimensional planar surface. Nevertheless, for small 
surface distances, one expects to recover the results for 
the two-dimensional scenario. The geometry furthermore 
represents a reasonable approximation to the tip of the 
needles of the Deslauriers et al. needle trap [M]. This 
is illustrated in Fig. [Tj where the three images from top 
to bottom represent the successive approximations made 
in going from the real experiment configuration to the 
geometry modeled here and in the following section. 

The Green's function for this system expressed in 
spherical coordinates defined in Fig. |8] is 



G(r, r') 



1 

47r 



1 



(28) 



Substituting into Eq. [T2j we obtain the following limiting 
forms for the geometric factor for an ion at distance r 



8 



from the center of a unit sphere (a = 1) 



A.(r) 
A.(r) 



, N 47r(T-2-l)* 

'o 

A 3(l+r^) 
, Af 47r(r2-l)* 



IP 
PP 

IP 

PP 



(29) 



Unhke the planar geometries, the geometric factor now 
has a non-trivial variation with distance. 

Although this geometry allows for a derivation of a 
in a closed-form, it is extremely helpful to consider the 
eigenfunction expansion of the geometric factor (Eq. 17 1. 



This approach lets us evaluate the geometric factor for 
intermediate patch sizes between the IP and PP limits 
(Eq. 18) from which we may gain intuition about how 



patch size affects the scaling of the power spectrum of 
electric field noise. Invoking this approach for a geometry 
that provides for a closed-form solution furthermore pro- 
vides for a consistency check between the two approaches. 

One could treat this problem rigorously by first solving 
Helmholtz's equation to obtain the eigenfunctions, but 
it is more expedient to start from the surface Green's 
function, which assumes the simple form: 

oo I i^i 

G.(r,r')=^ ^ ^Y,^{e',cb')Y;„,ie,cb). (30) 

/=0 m=-l 

Here the Yi„i are the spherical harmonics of degree / and 
order m normalized such that ^ ¥[^{0, 4')Yvm' <t))dVL = 
Su'Smm'- Substituting into Eq. [T2| 



A.(r) = 



A.(r) 



47r 

A 

N 



Q Yoo{e.,4.) 

r ^ 



A^ -r^oo v^/ 

N 2^1=0 Z^m=- 



d, 



„1+1 



IP 

PP 

IP 
PP 



(31) 



An advantage of the eigenfunction expansion of the ge- 
ometric factor in terms of spherical harmonics is that one 
may truncate higher order terms beyond some I — Iq > Q. 
Such an action roughly corresponds to evaluating the ge- 
ometric factor for patches of angular radii 9q ~ 2/Zo for 
large Iq ^ 0. Instead of sharp truncation, patch size de- 
pendance could also be introduced by attenuating higher 
order terms according to some distribution, resulting in 
different scaling behavior close to the surface. 

The scaling exponent a is evaluated, using both ap- 
proaches presented above, with respect to dimensionless 
distance D = d/a such that all dimensions are in units of 
electrode dimension, a. We make use of Eq.[T9]to extract 
a from the geometric factor and we perform the evalu- 
ation with respect to the surface of the unit sphere for 
both (r, 9) modes by setting Tq ^ d = f. Additionally, 
we exploit the spherical symmetry to set = </) = 0. Fig- 
ure [9] shows the scaling for the radial, r, mode as well as 
the transverse, 6, mode. 



The predictions in the limiting cases, given by Eq. 29 
are drawn in Fig. [9] as thick, solid lines, with IP and PP 
being the upper and lower lines, respectively. Similarly, 
the results of the approach based on the eigenfunction 



expansion, given in Eq. 31 are drawn as thin, dotted lines 
for patch sizes of 9^ ranging from 10^^ — 10^^ for the r 
and 9 modes as annotated in the figure. Note that for the 
9 mode, we have omitted the trivial case of truncating 
at Iq = 0, corresponding to a pure monopole with no 
component in the 9 direction. 

As expected, the predictions of the eigenfunction ex- 
pansion are all bounded by the solid lines representing the 
limiting cases of IP and PP. In the limit of D <^ 9(^, the 
results of the eigenfunction expansion converge towards 
the IP limit, since at infinitesimal distances, the finite 
patches appear infinite. The scale at which this occurs, 
however, is obviously patch size dependent. For example, 
if a = 1 cm, realistic patch sizes of 1 /im are obtained 
with 9q — lO"'* and the crossover between the IP and 
PP limits occur for surface distances of order 1 /im. This 
trend of decreasing a with decreasing distance scale D is 
in agreement with the infinite plane model of Dubessy et 
al. [21^. 

In the limit of _D 1, the scaling for all patch sizes 
converges toward 4 for the r mode and 6 for the 9 mode 
and become independent of the exact patch size, since 
at large distances, the ion is no longer sensitive to the 
structure of the surface. The result for the 9 mode is not 
obtainable in the infinite plane approximation. It is also 
worth noting that for _D ~ 10 - a scale commensurate 
with the Deslauriers et al. needle trap experiment [T3] 
- the scaling exponent for the r mode is consistent with 
the a = 3.5 ± 0.1 value obtained by that experiment. 



D. Spheroidal electrodes 

We now consider an ion suspended along the z axis 
outside electrodes of spheroidal geometries (Fig. 10). The 
spheroids approximate a wide range of shapes, and gen- 
eralize our discussion of spheres in spherical coordinates 
above. For example, a thin prolate spheroid in place of 
a sphere (Fig. [t]), would be a more reasonable approx- 
imation for the needle trap described in Ref. [M]. On 
the other hand, a finite planar trap, or disc, could be 
described as an extremely flat oblate ellipsoid. We may 
collectively model these geometries by solving Laplace's 
equation in spheroidal coordinates, where the geometry 
is effectively specified via a choice of either prolate or 
oblate coordinate system. 



We use the following definitions (c.f. Fig. 10) for pro- 
late spheroidal coordinates (1 <^ <oo,— 1 < rj < 1,0< 
(/> < 27r) 



X 



y = (^^ - 1)(1 - V'^) sin (/), 



z 



(32) 
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FIG. 9. Scaling exponent a as a function of dimensionless 
distance D = d/a above a sphere of radius a for both normal 
(r) and transverse {6) modes. The labeled solid lines are for 
the IP and PP limits, and are based on Eq. |29[ The shaded 
area between the IP and PP lines of each mode represents 
permissible values of a. The dotted lines correspond to in- 
termediate patches sizes eja = 10~^10"^10"^10"* eval- 
uated via Eq. 31 by truncating the expansion at la = 2/{0q). 
The dashed line of a = 4 is plotted for reference. 



and oblate spheroidal coordmates (0<^<oo,— 1 < t] < 
1,0 < < 27r) 



X — a\J (1 + ^^)(1 — ?7^) COS0, 
y = aV(l + C^)(l-?7^)sin(/>, 
z = a^r]. 



(33) 



In this notation, the r] and ^ coordinates are the 
spheroidal analogues of the spherical coordinates cos 61 
and r respectively. For example, in prolate spheroidal 
coordinates, surfaces of constant ^ = form elongated 
ellipsoids. The degenerate limit ^ 1 is a straight 
line between z = ±a, while the other limit ~^ oo is 
a sphere. In oblate spheroidal coordinates, surfaces of 
constant ^ = Co form flattened ellipsoids, and the limits 
Co — and Co ~^ oo correspond to a disc and a sphere, re- 
spectively. For simplicity, we consider the unit spheroid, 
meaning that all dimensions are in units of a. 

Unlike in the analysis for planes and spheres, the 
Green's function for these geometries cannot be written 
in closed form. However, Laplace's equation is separable 
in the spheroidal coordinates, and the general solution 



may be expressed as a series of Legendre polynomials, 
Pimix) and Qim{x), of the first and second kind, respec- 
tively, defined with normalization 



Plm{x)Pvm{x)dx = 5w 



(34) 



Hence, similar to the eigenfunction expansion of the geo- 
metric factor, used in our treatment of spheres, we start 
from the surface Green's function for prolate and oblate 
coordinates. We then examine the scaling of the geomet- 
ric factor Afc(r) in Eg. [I2| along the z axis, where -q = 1. 
a) z| b) zi 




FIG. 10. (a) Prolate spheroidal coordinates defined in Eq. 32 
(b) oblate spheroidal coordinates defined in Eq. |33] 



Both at very small and at very large distances from the 
electrodes, we intuitively expect scaling laws similar to 
those found for the sphere; however, at intermediate dis- 
tances, the shape of the electrodes is expected to become 
important. 

Given an ellipsoid with surface Coj the prolate spheroid 
surface Green's function is 



G.(r,r') = 



and the oblate spheroid surface Green's function is 



(35) 



G,(r,r') 



2vrv/(r+c!)(cf+?^ 



, (36) 



where m = —I, — ^ + 1, • • • ,1 in the sum. By substitut- 
ing Eq. [35] and [36] into Eq. [12] the geometric factors for 
prolate and oblate spheroids become 
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FIG. 11. (a) Needle-shaped electrode in prolate spheroidal 
coordinates for a surface ^ > 1. (b) Disc-shaped electrode in 
oblate spheroidal coordinates for a surface ^ = 0. The ion is 
represented by the black dot. 



where the sum over m is for |m| < min(Z,^') and where 



-dry, 



and 



QUO Pimivy""^ 

Qimi£,o) \/27r 



(39) 



(40) 



We note that these results reduce to those of the sphere 
in the limit — J' oo as expected, since limj_j.oo QimiO 

We evaluate the geometric factor numerically for the ^ 
and r] modes along the z axis. In this direction, all terms 
with m ^ (|m| 1) vanish for the ^ (r;) mode. The 
evaluation is complicated by the cross terms c^,^ as the 
system is not invariant with respect to the t] coordinate. 
However, not all terms are significant as c^,^ is observed 
to go as ~ e"l'"'V«o-i/«o foj. j^rge |/ - 

The geometric factor, Afc(r), is computed numerically 
for both needle and disc geometries by taking appropriate 
limits for prolate and oblate spheroids, respectively. We 
perform the sum over 1,1' by truncating a,t I + I' — 21q. 
Similar to the case of the sphere, we define truncation at 
^0 to correspond a patch size of 9q — 2/Iq. The lowest 
order non-vanishing term for the (77) mode occurs at 
Iq = {Iq — 1). In both cases, scaling is evaluated with 
respect to distance d along the z axis, as shown in Fig. [TT] 



1. Prolate spheroid: needle 

The general features of the geometric factor for a 
prolate spheroid are illustrated through an example ge- 
ometry with an aspect ratio, half-length/radius — 100. 
The high aspect ratio mimics that of needle-shaped elec- 
trodes as illustrated in Fig. [t] (a) and (b). Specifically, 
we model the coordinate surface ^0 = 100/(3\/TTTT) ss 
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5 
4 

ii 3 
2 
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FIG. 12. Scaling exponent a as a function of dimensionless 



distance D 



d/rpr, 



above a needle of radius r 



prolate 



0.0100005 X a and half-length 100 x rpioiato for both normal 
(^) and transverse (rj) modes. The labeled solid lines are for 
the IP limits. The dotted lines correspond to intermediate 
patches sizes 8q = 0.4, 0.04 evaluated via Eq.[37]by truncating 
the expression at I + I' — 2lo = 4/(Sj). The dashed line of 
a = 4 is plotted for reference. 



1.00005. This represents a spheroid with a half-length of 
^oflproiatc ~ al. 00005 and a radius, rpi-oiatc ~ aO. 0100005, 
at the widest point. All dimensions are modeled in units 
of a and we set Tq = ^o-^, d = z, 77 = 1 in the geo- 
metric factor of Eq. [12] With this choice of geometry 
^ = z, and d refers to distance above the top of the 
needle [c.f. Fig. [Tl] (a)]. The dependence of the scaling 
exponent, a, with respect to the dimensionless distance, 
D = (i/(rproiatc), for the ^ and 77 modes are plotted in 
Fig. [121 

The qualitative behavior of the scaling exponent for 
a prolate needle largely imitates that of the sphere: In 
the limit D ^ 1, the ^ mode converges to a = 4 and 
the 7] mode to a = 6, while in the limit oi D ^ 1, both 
modes exhibit convergence toward a — 0. Compared to 
the results for the sphere, however, additional features 
have arisen. The brief plateau in a that occurs for inter- 
mediate 0.1 < Z? < 100 is a geometric effect due to the 
elongated shape of the needle. For lower aspect ratios 
we find that the region of this plateau narrows and that 
convergence toward a = becomes evident already at 
larger values of D. We note that the bump in a, for the 
^ mode around D ^ 1, is also present in the exact ana- 
lytic solution for spheres, and is not a numerical artifact. 
Such non-nionotonic behavior is not entirely surprising 
as this region represents the cross-over between the two 
distinct regimes of infinite and infinitesimal ion-surface 
distance. 

Similarly to the results of the sphere, we find, for the 
^ mode in the region explored by the Deslauriers et al. 
needle experiment [14 , that the possible values of a are 
bounded to an interval below 4 but consistent with the 
a = 3.5 ± 0.1 reported by their experiment. 
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FIG. 13. Scaling exponent a as a function of dimensionless 
distance D — d/a above a disc of radius a for both normal 
(^) and transverse (ri) modes. The labeled solid lines are for 
the IP limits. The dotted lines correspond to intermediate 
patches sizes 9q = 0.4, 0.04 evaluated via Eq.[38]by truncating 
the expression at I + I' = 2lo — 4/(6^). The dashed line of 
Of = 4 is plotted for reference. 



2. Oblate spheroid: disc 



The most extreme oblate spheroid is a disc with a ra- 
dius of — a and zero thickness. This geometry may 
be used to model the surface electrode ion trap 23 25J , 
currently receiving much attention due to its potential 
in quantum information science |35j - a field that expe- 
riences particular sensitivity to electric field noise |15j . 
Considering this geometry furthermore extends the work 
of Dubessy et al. ^21j to finite planes. 

For our numerical evaluation, we set rp = 0, d = z, 
and rj — 1 in the geometric factor of Eq. [12] With this 
choice of geometry, d refers to distance above the disc 
origin [c.f. Fig. [Tl] (b)]. A complication occurs when 
evaluating a for the disc: As in the case of the hole trap, 
where the thin edge caused a logarithmic divergence in a, 
a similar effect arises here. This is resolved by restricting 
the region of integration over r] in the coefficient to 
(— 1, —6), ((5, 1), where < (5 <C 1 to avoid the edge. Here 
we have used 5 = 0.1. 

The scaling exponent, a, for the ^ and 77 modes of 
the disc are plotted in Fig. [TSj The qualitative behavior 
is similar to the prolate needle geometry and, in some 
limits, to the infinite plane: In the limit of Z? = d/a <C 1, 
e.g., convergence toward a = is observed, as expected 
in this limit where all surfaces appear infinite regardless 
of their exact geometry. In the limit of D 3> 1 the ^ 
mode converges to a = 4, while the 77 mode converges 
toward a = 6, once again signifying a strong departure 
from the results of the infinite plane model. 



IV. CONCLUSIONS 



We have presented an analytic model for electric field 
noise from fluctuating patch potentials starting from 
Laplace's equation and its Green's function solution. Be- 
yond geometries for which the Green's function is readily 
obtained, our model uses an eigenfunction expansion and 
we employ this to analyze geometries of relevance to cur- 
rent ion trapping technology. 

At distance scales that are relevant to typical ion traps, 
prior works have collectively established a scaling for elec- 
tric field noise with surface distance, d, of d^", with 
a — 4. While our model is in agreement with those 
models in a number of scenarios, for certain parameter 
regimes we observe a strong dependence on the finite ge- 
ometry as well as on the patch size, leading to a departure 
from the a = 4 scaling. Moreover, we consider the effect 
of electric field noise on motional modes of a trapped ion 
that are both normal and transverse relative to the elec- 
trode surface studied. Significant differences between the 
two orthogonal modes are predicted for the more extreme 
geometries such as e.g. needle shaped electrodes. 

It would be of interest to confirm these predictions 
experimentally, and a suitable geometry could be either 
the needle trap of Deslauriers et al. [13] or the recently 
developed stylus trap, which has also been proposed as 
a highly sensitive electric field probe [35]. Combined 
with techniques for varying the ion-surface distance, d, 
in-situ [371 EH] 7 such systems could be used for detailed 
tests of our model. 

Although the geometries we have considered are mostly 
finite, not all are representative of real ion traps. In 
particular, the needle trap of Deslauriers et al. [14j was 
approximated by a single needle rather than two as in 
the actual experiment. To what extent this influences 
the scaling is difficult to gauge intuitively, but it could 
potentially be investigated by modeling two spheres ex- 
actly in bispherical coordinates where Laplace's equation 
is separable. 

More generally, the geometries considered in this work 
have been generic shapes for which solutions to Laplace's 
equation can be written out either through the appropri- 
ate Green's function or via an eigenfunction expansion. 
Real ion traps often differ from such geometries; however, 
if the electrostatic Green's function can be obtained nu- 
merically, the approach presented in this work could be 
extended to arbitrary geometries. In this respect we em- 
phasize the method of the eigenfunction expansion. As 
it in principle allows for the study of arbitrary patch cor- 
relation functions, it may find use in modeling of exper- 
imental scenarios relevant, for example, to Casimir force 
measurements and nanoscale surface probes, where the 
characteristic scale of the probe-surface distance becomes 
commensurate with the effective patch size. 
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